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1 Introduction 



The indirect CP violation, present in the neutral Kaon meson pair (K° — k") os- 
cillations pj, is a confirmed experimental fact. The neutral Kaons K° = {sd) and 

K = (sd) are flavour eigenstates which, in the Standard Model (SM), can mix to 

lowest order through a box diagram due to the weak interactions. The weak hamil- 

tonian describing the time evolution of the kaon pair system can be expressed via 

a 2 X 2 non-hermitian matrix, H = M — iV /2. The dispersive and absorptive 

parts, M and F respectively, are hermitian matrices. Defining the CP eigenstates 

l-^it) = (1-^*^) i \K ))/a/2, it turns out that the two hamiltonian eigenstates, called 

i^— short and i^— long, are I-R's'/l) = / ^ -o (l-^±) + ^I-^t))- The paramater e 

V 1 + e'' 

represents the mixing between the two CP eigenstates; it depends on the phase con- 
ventions of the eigenstates and represents the theoretical measure of the indirect CP 
violation. On the other hand the experimental measurement involves a phase inde- 
pendent quantity which is 

_ A{Kl (7r7r),=o)) 
^ A{Ks (vr7r),=o)) ^ ^ 

It turns out that (Ref. p]) ex = e + i^, where ^ = j^^j^^ and Aq is the isospin zero 
amplitude for the K" decay to either (vr^vr^) or (tt^tt'^). It has been shown ( [3],[1], 
[5J; see also [6]) that, under reasonable estimates for the long-distance contribution to 
both the dispersive and the absorptive parts of the hamiltonian, ex is approximately 
given 

e^'^^ImMia 

In Eq. (E]) the phase 0^ = 43.51(5)° and the mass difference between Kl and Kl, 
Am-K = 3.483(6) x 10~^^ MeV, have been determined experimentally with high pre- 
cision [7]. The factor = 0.94(2) incorporates in an approximate way the effect of 
the long distance contributions mentioned above. 

In the SM the calculation of Im M12 is performed in the Operator Product Expan- 
sion (OPE) approach. Once the heavy degrees of freedom, including the charm quark, 
are integrated out the CP violation contribution in the neutral K-meson mixing is 
described in terms of the K — matrix element of a local AS" = 2 four-fermion 
operator: 

Mu = ^((^l^f/r 1^°))* = ^C^iMK'\0^'='ifi)\K')r (3) 

The Wilson coefficient Cy^(/x) contains all the short distance effects which are cal- 
culated in perturbation theory. The matrix element of the four-fermion operator 
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O ^ = Ci = {sd)v-A{sd)v-A = s'yi_i{l — 75)6? S7^(l — 75)0? is written in the form 

(K'm^^)\K') = ^mylBM (4) 

where the parameter Sk is the amount by which it differs from its Vacuum Saturation 
Approximation (VSA) value and mx and /k are the mass and the decay constant of 
the Kaon respectively. Therefore the bag parameter Sk is the measure of the non- 
perturbative QCD contribution in the hadronic matrix element of the K-meson pair 
oscillation. The Renormalisation Group Invariant (RGI) quantity is defined by the 
equation 

1-2/9 r a'f\fi) 



tt(3)(^) i + r^j3 5j,(^) (5) 



where (/x) is the running coupling constant for Nj = 3 flavours and J3 has been 
computed up to NLO order |8j. Eqs. (E])-© yield 
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lex] = KeC.BK iKfepA^r] - r]iSo{xc)il - + '^■^^oixt)\Vcb\'^X\l - p) + r]sSo{xc, Xt 

(6) 

where = '^^/^T^^^ ^^"^ — iKisI- 'S'o are the Inami-Lim functions depending 
on Xc^t = "^ctf-^w ^^"^ Sive the charm, top and charm-top contributions to the box 
diagram; ?7i=i,2,3 contain the corresponding short distance QCD contributions to NLO 
ordeiE ([9]-[n]'). The experimental value of je^l = 2.228(11) X 10 ^ is known with 
high precision [7]. Therefore, the hyperbola defined by Eq. ([6]) in the (p, rj) plane can 
be used to (over-) constrain the upper vertex of the unitarity triangle with a precision 
that depends on the quality of the estimates of \Vus\, \Vcb\ and -Bk- While the nu- 
merical value of the first is known to very good precision (see Ref. [I3]), the estimate 
of (inclusive) \Vcb\ (Ref. [T^]) is still given with un uncertainty of ~ 2% which gets 
amplified four times in Eq. (|6]) because it enters to the fourth power. can be 
calculated on the lattice from first principles. Its estimate, until 2008, used to repre- 
sent the largest source of uncertainty in Eq. (E]) but now, thanks to recent precision 
lattice calculations, it is known with a total (statistical and systematic) uncertainty 
of less than 4%. This is the consequence of employing unquenched simulations with 
Nf = 2 and Nf = 2 -|- 1 dynamical quarks with high statistics, controlling with better 
accuracy both the discretization errors and the extrapolation to the physical point 
and using non-perturbative methods for the renormalisation of the operators on the 
lattice. 



"See Ref. [12] on a recent calculation for 773 to NNLO order. 



2 



2 Bk calculation on the lattice 



2.1 General considerations 

The evaluation of 5k on the lattice requires the computation of a three point corre- 
lation function with the insertion of two pseudoscalar meson interpolating fields at 
two time slices, say, ti and t2 and the insertion of the AS = 2 four-fermion operator 
considered at any time slice to with to G [ti, ^2]- In order to carry out the computation 
of Eq. dl]), the three point correlation function needs to be divided by two two-point 
correlation functions each of which involves a pseudoscalar meson and an axial cur- 
rent interpolating fields; then one takes the asymptotic limit i.e. ti ^ to ^ ^2 in 
order to obtain the value of the matrix element between the lowest-lying states. In 
summary, one has 

(8/3)(P(ti)tA(to))(A(to)tP(t2)) i8/3)(K'\Ao\0){0\Al\KO) 

where V(t) and A(t) denote the pseudoscalar density and axial current interpolating 
fields respectively. Note also that (OlAol-ft'^) = /k^k- In practice, due to parity 
conservation in QCD, one needs to evaluate only the matrix element of the parity- 
even part, Ovv+AA = sj^ds'-y^d + 57^75^57^756?, of the O^^^"^ operator. 

Naturally, a renormalization step is necessary to get finite results in the continuum 
limit (a — )■ 0). For this purpose the renormalization constant (RC) of the operator O, 
Zo, (or a full matrix of RCs in the general case where more than one operators get 
mixed in the renormalisation procedure) has to be computed at a certain scale /i (the 
same as that of the Wilson coefficient). Lattice perturbation theory can be used. The 
result, however, suffers from more or less insufficiently well estimated systematic errors 
due to the truncation of the perturbative expansion. This is especially true for not 
so fine lattice spacings. In any case non-perturbative methods offer a more accurate 
way in computing the RC of operators. By imposing renormalization conditions on 
correlation functions with operator insertions and performing calculations of matrix 
elements on the lattice, one has the opportunity to embody in the calculation both 
higher order contibutions and non-perturbative effects. In this way the remaining 
sources of systematic error are only due to the use of perturbation theory in passing 
from the lattice scheme (at some scale fi) to a continuum renormalisation scheme. 

Lattice calculations are performed at a number of various values of lattice spacing 
and suffer from discretization errors which can be eliminated when the extrapolation 
to the continuum limit (c.l.) is achieved. This source of systematic error can be well 
kept under control if the values of the lattice spacing are small, typically less than 
about 0.1 fm. Then simulations for at least three values of the lattice spacing are 
typically necessary to be able to perform the extrapolation to the continuum limit. 
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In most of the cases discretization error of the various physical quantities measured 
on the lattice are of at most O(a^). 

It is interesting to have a look at the "world-average" values of Bn which cover 
a period of the past fifteen years as they reported at the Lattice Conferences. The 
RGI results (computed with the same running using Nf = 3 flavours, see Ref. [T5] ) 
are shown in Fig. [H 



' LAT1996 LAT2000 LAT2005 LAT2008 LAT2009 CKM2010 

Figure 1: i?K "world averages" over the last fifteen years reported at the Lattice 
Conferences. The world average result at CKM2010 is also included. 

It is impressive to notice the great improvement in the quality of the results 
where total uncertainty has decreased by almost five times in the last fifteen years. It 
should also be noticed that even until 2005 the Bk "world average" result was based 
on estimates computed in the quenched approximation with the systematic error due 
to quenching being in practice unknown. However, the "world average" value of the 
last two years is based only on results produced with unquenched simulations. 

2.2 Recent lattice Bk computations 

Different lattice groups compute Bk using a variety of lattice regularizations which 
conserve (most of) the chiral symmetry and provide well controlled and small O(a^) 
discretization errors. Besides that, simulations have been performed at a well tuned 
physical value for the strange quark mass while the lightest simulated quark mass 
value, even if not yet at its physical value, has smaller values than in any preceding 
calculation of -Bk- Four collaborations, performing unquenched simulations, have 
computed the continuum limit value for Bk- Three of them have simulated Nf = 2 + 1 
dynamical quarks, while the fourth one has used Nf = 2 dynamical quark simulations. 

Table [T] gives a summary of the simulation details for the computation by each 
collaboration; Table [2] contains B^ results expressed in the MS scheme at 2 GeV as 
well as in the RGI definition (see Refs. [IE]- [19]); the first error is statistical while 
the second is the systematic one; the total error has to be calculated as the sum in 
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Collaboration 


Nf 


fermion discr. 
val. / sea 


a [fm] 


r , 1 

L J min 

in [MeV] 


r , 1 

\mt'l\ . 

L J nun 


Renorni/tion 


ALV [E] 


2+1 


DW/Asqtad 


0.12, 0.09 


240/370 


3.5 


Non-Pert. 


BSW [IT] 


2+1 


HYP-stag/Asqtad 


0.12, 0.09, 0.06 


240/300 


2.5 


Pert. 


RBC-UKQCD [IB] 


2+1 


DW/DW 


0.11, 0.09 


220/290 


3.1 


Non-Pert. 


ETMC PJ 


2 


OS/TM 


0.10, 0.09, 0.07 


280/280 


3.3 


Non-Pert. 



Table 1: Simulation details of unquenched simulations of Bk computation with Nf = 
2 and Nf = 2 + 1 dynamical quarks. 



quadrature of the statistical and the systematic ones. In the following paragraphs a 
short description of each computation is presented. 

Aubin, Laiho and Van de Water (ALV) (see Ref. (TB]) have presented a mixed 
action calculation performed with Domain Wall (DW) fermion regularization in the 
valence quark sector on Nf = 2 + 1 dynamical Asqt ad- improved staggered quarks 
(MILC). DW fermions in the valence sector offer the possibility to use the MILC 
dynamical configurations without having to bother about taste mixing. Since they 
have a residual chiral breaking (due to the finite size of the fifth lattice dimension), 
the renormalisation of the parity even four-fermion operator suffers from mixing with 
operators of wrong chirality through coefficients whose values are expected to be 
suppressed at order 0((amres)^)- The implementation of the non-perturbative renor- 
malisation method (RI-MOM) is straightforward for DW fermions. They have done 
simulations at two values of lattice spacing, namely a = 0.12 and 0.09 fm. They 
extrapolate to the physical point using fit functions based on Mixed action xPT with 
the addition of NNLO analytical terms. Their result is given in Table [2] ; a large 
part of the systematic error quoted is due to the uncertainty in the renormalisation 
constant computation. 

Bae et al. (BSW|l], see Ref. [T^) have also completed a mixed action computation 
where they have used HYP-smeared staggered valence quarks on MILC 2-1-1 dynam- 
ical quark configurations. The choice for this particular valence quark regularization 
is justified by its property of dispalying rather small taste breaking effects and a com- 
putationally cheap implementation. They extrapolate to the physical point trying 
fit functions based on SU(3) and SU(2) Staggered-xPT. The two fitting procedures 
give compatible results with the latter being more straightforward in fitting the data. 
Simulations have been carried out at three values of lattice spacing, a = 0.12, 0.09 
and 0.06 fm. The dominant source of uncertainty is due to the perturbative approach 
they use to renormalize the four-fermion operator. In fact, the total error is almost 
6% the larger part of which is due to RC's uncertainty. 

RBC-UKQCD Collaboration (Ref. [IB]) used a iV/ = 2 + 1 DW action both for 

'It is an acronym for Brookheaven, Seoul and Washington groups. 
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the sea and the valence quarks at two values of the lattice spacing, a = 0.11 and 
0.09 fm. The renormalization of the four-fermion operator is performed employing 
RI-MOM methods where non-exceptional momentum renormalization conditions and 
twisted boundary conditions have been applied. They achieve further suppression of 
the residual wrong chirality mixing and a matching to MS at a higher scale value (3 
GeV) which allows for a a better control of the perturbative systematic uncertainty. 
The uncertainty to the final value due to the renormalisation procedure is estimated 
2%. The physical point is reached via a combined continuum and SU(2) PQxPT 
fitting formula. The total error in the final value is 3.6%, while the pure stastistical 
one is about 0.9%. 

ETM Collaboration (Ref. [19]) has adopted a mixed action set-up to perform the 
computation with a Wilson fermion regularization. Actually, since Wilson fermions 
explicitly break the chiral symmetry, the relevant computation done with the use 
of plain Wilson quarks suffers from 0(a) discretization errors and yields wrong chi- 
rality mixing in the process of renornalizing the lattice Ow+aa operator. It is known 
that both problems can inject serious systematic uncertainties in the computation. 
However, it has been shown that they can be tackled simultaneously by employing 
Osterwalder- Seller (OS) valence on Twisted Mass (TM) sea fermions, both tuned at 
maximal twist (see Ref. [2Q]). ETMC has carried out the computation working with 
Nf = 2 dynamical quark simulations, so the strange quark is still quenched. Unitarity 
violations due to the mixed action set-up, that are expected to be O(a^) effects, have 
been shown to be well under control in the continuum limit. The extrapolation to 
the physical light quark mass is performed using a fit function based on SU(2)-xPT 
while three values of the lattice spacing, a = 0.10, 0.09 and 0.07 fm, have been used 
to extrapolate to the continuum limit. The renormalization of the four-fermion oper- 
ator has been carried out in a non-perturbative way using the RI-MOM method. The 
mixing with operators of wrong chirality has been shown to be numerically negligible. 
The total error is 4%. 

The i?K results published in the last five years and computed on quenched {Nf = 0) 
and unquenched lattices with Nf = 2 and Nf = 2 + 1 dynamical quark simulations 
are collected in Fig. [2J The label (c./.) is attached to the 5k values obtained after the 
continuum limit has been taken. Vertical lines indicate the average result with 
its error, based on Nf = 2 + 1 simulations, which is 

= 0.736(05)(26) @ Nf = 2 + 1 (8) 

This value has been obtained by taking the average of the Nf = 2 + 1 results weighted 
with the quoted statistical error. The systematic uncertainty is the smallest error 
quoted in Table [2l The total uncertainty is about 3.6%. The Bk average for the 
Nf = 2 case obviously coincides with the result provided by ETMC, since it is the 
only continuum limit result available. One finds 

5k = 0.729(25)(16) @ Nf = 2 (9) 
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Collaboration 




AT 


-Dj^ l^Z Vore V ) 


D 

-DK 


A T \^ 
AL/V 




Z + i 


U.OZ ( [Uv ) [A\J ) 


n 70/1 ^ns^ ^OQ^ 
U. ( /4(^Uo J (^/o J 


BSW 




2+1 


0.529(09)(32) 


0.724(12)(43) 


RBC-UKQCD 


[T8l 


2+1 


0.549(05)(26) 


0.749(07) (26) 


ETMC 


m 


2 


0.517(18)(11) 


0.729(25)(17) 



Table 2: Results in the continuum limit for the -BK-parameter expressed in the MS 
scheme and in the RGI definition from unquenched simulations with Nf = 2 and 
2 + 1 dynamical quarks. The first error is quoted as statistical while the second is the 
estimate of the systematic uncertainties. 

Results from Eqs. ([8]) and Qj indicate that the systematic uncertainty due to the 
quenching of the strange quark has an impact which is smaller than other system- 
atic errors. Moreover, comparing the result of Eq. (|8]) with the most precise of the 
quenched results (see Ref. |25j), one can infer that the systematic uncertainty due to 
the quenching is small (actually ~ 6%). 

mixing beyond the SM 

There are various models for the physics Beyond the SM (BSM) which lead to other 
possible AS = 2 processes at one loop. In this case the computation of the relevant 
matrix elements of the effective hamiltonian in combination with the experimental 
value of ex would offer the chance of obtaining constraints on the parameters of the 
model (like for instance estimates on the off-diagonal terms of the squark matrix in 
supersymmetric models [22]) which enter explicitly in the Wilson coefficients. 
The general form for the AS* = 2 effective Hamiltonian in BSM models reads 

5 3 

nti=' = $^a(/i)0. + J2c^{fi)o^ (10) 

i=l i=l 



where 








01 = 








O2 = 


r(l-75)c?1[^''(l-75)ci1, 


03 = 


r(i-75)rfi[^''(i-75)rf" 


04 = 


ini - j,)d^][s\i + j,)d% 


05 = 


[r{i-^,)dV{i + ^,)d'^[ 


01 = 


[r^^{i + ^,)d'^][s\{i + ^,)d% 






02 = 


+ j,)d^][s\l + J,)d% 


03 = 


[r{l + ^,)d'][s\l + ^,)d^] 
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.......................... 


[231 Y Aoki et al f2006) h 


. 1 nw (c I ) 


[21] ALPHA (2007) i . 


, 1 TM (cJ.) 


[25] CP-PACS (2008) 


H DW {c.l.) 


[22] ALPHA (2009) M 


H OS-TM {c.l.) 


Nf = 2 




[261 Y Aoki et al f2005) han 


DW 


[24] JLQCD (2008) h- 


OV 


[19] ETMC (2010) hi 


H OS/TM (c./.) 


Nf = 2 + 1 




[16] ALV (2009) 


H DW/Asqtad (cl.) 


[17] BSW (2010) ^ 


H HYP-stag/Asqtad (c./.) 


[18] RBC-UKQCD (2010) h 


^ DW (cl.) 



0.6 0.7 0.8 0.9 



Bk 

Figure 2: 5k results from Nf = (quenched) and unquenched simulations with Nf 
= 2 and Nf = 2+1 dynamical quarks. Computations where the continuum limit has 
been taken are indicated with "(c./.)". Vertical lines refer to the average result of 
Eq. ([H]). In the plot, results from pioneer works with A^^ = 2 [27] and Nf = 2+1 [2S] 
unquenched simulations have not been included since their quoted total uncertainty 
is now considered too large compared to today standards. 



We have seen that in the SM case, Eq. ([3]), only the operator Oi contributes. 
The parity-even parts of the operators Oi coincide with those of the operators Oi. 
Therefore, due to parity conservation in the strong interactions only the parity-even 
contribution of the operators Oi need to be calculated. A basis of the parity-even 
operators is 

= is^5d){s^5d), O'' = {sd){sd), 

O^^ = isa^,d){sa^,d) (12) 

in terms of which and after using a Fierz transformation one obtains the relations 

O, = (O^^ + O^-^), 02 = (0^^ + 0^^), 03 = -i(0^^ + 0^^-0^^), 
O, = (O^^-O^^), = -1(0^^-0^^) (13) 



The B-parameters for the operators of Eq. f[T3|) are defined as 
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where Cj = {—5/3, 1/3, 2, 2/3}, z = 2, . . . , 5. The matrix element of the operator Oi 
vanishes in the chiral hmit while the matrix element of the operators Oi, i = 2, . . . ,5 
get a non-zero value in the chiral limit. From the above equations it can be seen 
that the calculation of the Bi parameters for i = 2, ... ,5 involves the calculation 
of the quark mass at a common renormalization scale fi. In order to avoid any 
extra systematic uncertainties in the computation of the matrix elements due to the 
quark mass evaluation, an alternative calculation has been proposed which consists 
in calculating directly appropriate ratios of the Oi, i = 2, . . . ,5 four-fermion matrix 
elements with the Oi onS (pQ], |3T]). 

There are a few quenched lattice calculations performed with the use of tree-level 
improved Wilson fermions [3D] , overlap fermions [21] and DW fermions [52] ■ However 
only the first two computations have been carried out using two values of the lattice 
spacing allowing for the possibility of estimating discretisation effects. A preliminary 
study referring to the bare matrix elements computation using Nf = 2 + 1 dynamical 
DW fermions was also presented in Ref . [33] . Recently ETMC has presented a Nf = 2 
unquenched calculation using the mixed action set-up described above (Section 2.2). 
Both Bi-hag parameters and matrix elements ratios have been calculated at three 
values of the lattice spacing, a = 0.1, 0.09 and 0.07 fm from which reliable continuum 
limit estimates can be obtained. The results presented in Ref. [31] are preliminary; 
however, it might be useful to proceed at a first comparison between the quenched 
and the unquenched estimates for the bag-parameters Bi, i = 2, . . . ,5 calculated in 
the MS scheme at 2 GeV. This comparison is provided in Fig. [3l 



4 Conclusions 

Bk computation on the lattice has already entered the era of precision measurements 
with fully unquenched simulations. A number of lattice collaborations using various 
lattice regularizations provide results that are in nice agreement among themselves. 
The total estimated uncertainty is less than 4%; as a consequence the lattice cal- 
culations are not responsible any more for the main source of uncertainty in the eK 
equation (Eq. |6]). The comparison of B^ results produced with Nf = 2 + 1 and 

■l-Note that in forming these ratios one should take care of the fact that the matrix element in 
the denominator vanishes in the chiral limit. 
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Figure 3: Comparison of continuum limit Bi 



2, . . . , 5) results (in MS scheme 



at 2 GeV) obtained using Nf = {] (quenched) ([30], [3T]) and unquenched Nj = 2 
simulations [34j. 



Nf = 2 unquenched simulations indicate that the systematic error due to the quench- 
ing of the strange quark is smaller than other systematic uncertianties. The variety of 
fermion regularisations, today in use, offer a very good control over the contamination 
of wrong chirality operator mixing in the renormalisation of the AS* = 2 four-fermion 
operator. Multiplicative renormalization is guaranteed with discretization errors of 
O(a^). Many collaborations use non-perturbative methods for the operator renormal- 
ization leading to a good control over an important source of systematic uncertainty. 
However the operator renormalisation still representes about half of the systematic 
uncertainties. Another source of systematic error comes from the fit extrapolation to 
the physical point; both the choice of the fitting function and the relatively large value 
of the lightest quark mass actually used in the simulations for 5k inevitably inject 
a non- negligible systematic uncertainty into the final result. Of course simulations 
at finer lattice spacing in the near future would make smaller the systematics due to 
lattice artifacts making more reliable the continuum limit extrapolation itself as well 
as the evaluation of renormalization constants. 

Recently a Nf = 2 unquenched computation of the matrix elements of the full 
AS = 2 operator basis emerging from box diagrams in models beyond the Standard 
Model has been presented by the ETM Collaboration using three values of the lattice 
spacing. The continuum limit extrapolated results offer the opportunity of providing 
constraints on the input parameters of a class of models BSM. 
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